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FOREWORD 

Tlic Aeronomy Laboratory is organized into working groups according to 
the various areas of research being carried out. This arrangement is intended 
to encourage better communication among graduate students and staff members, 
and to stimulate ideas for new lines of scientific inquiry. Individuals may 
participate in as many groups as their interests dictate. 

This Progress Report for the period October 1, 1980 - March 31, 1981 con- 
tains summaries of work performed by each working group. These summaries 
have been written as a part of the groups' activities. 

On the following page the various groups and members are listed in the 
order that their reports appear in this volume. Support for the research 
of each group is given in parentheses following the name. 
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1. IONOSPHERIC 1HEORY 

This project involves theoretl''al and experimental studies of the effects 
of atmospheric dynamics on the ionosphere and neutral chemistry of the strato- 
sphere, mesosj^ere and thermosphere. Tlie research is supported in part by 
the National Aeronautics and Space Administration utdor grants 14-005-181 and 
NSG 5213, and in part by the National Science Foundation under grant A’I>1 79- 
09982. 

1.1 Planetary-Wave Dynamioe 

An investigation i,s continuing into the effects of vertically propagating 
planetary waves on NO densities in the 70 to 120 km altitude region. A time- 
dependent model has been developed to calculate the transport of NO due to 
planetary-wave- induced winds. These transports seem to be sensitive to the 
plnnetary-wave forcing used at the lower boundary, as well as to the structure 
of the zonal mean basic state. 

The results shown are from a calculation of steady-state NO transport 
by stationary planetax'y waves. Fluxes and transports of NO are calculated 
for the months of .January and February, 1974, using a separate basic state 
and '»oundary forcing for eacli month. Tixe boundary forcings were obtained 
by performing a spatial Fourier decomposition on .4 mb NMC geopotential height 
data for the weeks of January 2 through February 27, 1974. The zonal mean 
wind profile for the months of January and February were also inferred from 
NMC data. The month of February was interesting in that a major stratospheric 
warming occurred near the end of that month. The large planetary wave 
amplitude and reduction in zonal mean wind velocities due to this warming 
are reflected in the February average boundary forcing and zonal mean state. 

Figure 1 . 1 show s the results o f a cal culation of zonally averaged 
horizontal (i»* [NO) ' ) and vertical (u?' [NO] ' ) NO flux at 105 km. An upward 
arrow represents a northward upward motion. These fluxes are due to plane- 
tary wave number 1 for the month of January and February. It can be seen 
that while the horizontal flux is to the south in January, the directi'’n 
reverses and the flow is to the north in February. The vertical flux also 
undergoes this same reversal, with upward flux in January and downward flux 
in February. 

By calculating vertical and horizontal transports, the effect of these 
fiu,’.es on the zonally averxiged NO concentration at 105 km can be determined. 
The zonally averaged continuity equation for NO in the absence of a mean 
meridional flr>w can be expressed as: 

- 1 ) 

The first and second terms on the right hand side represent the change in the 
zonally averaged nitric oxide concentration due to north-south and vertical 
transport of nitric oxide by planetary waves. 
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Figure 1.1 


Zonal ly averaged horizontal and vertical NO flux 
for wave number 1 at 105 km. 
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Figure 1.2 showi*. the change In [NO] due to north-south transport by wave 
number 1 for the months of January and February. In January, the concentration 
is decreasing in the polar region and increasing in midlatitudes, due to a 
southward flux of NO. In February, however, the situation is reversed, with 
the 2 ,onally averaged concentration increasing in the polar regions and de- 
creasing in midlatitudes duo to a northward flux of NO. 

Figure 1.3 shows the ch.mge in NO due to vertical transport by wave 
niunber 1. Note that the changes in Figure 1.3 are an order of magnitude 
smaller than those in Figure 1.2. The reversal from month to month ‘‘ecn in 
Figure 1.2 is also apparent in Figure 1.3. In January, the concentration 
of NO at 105 km is increasing due to a flux of nitric oxide from below. 

In February, however, the concentration at 105 is decreasi.ig, due to a flu; 
of NO downward out of the region. 

A more detailed investigation into the area of NO transport is current!/ 
underway. A time -depen dent model will be used to calculate the transport 
of NO by both planetary waves and the mean meridional circulation. The 
effects of changes in zonal ly averaged production and loss rates of NO will 
also be considered. 


CHANGE IN [no] WITH T;mE IO^(l/sl(l/cm*) 


4 




CHANGE IN [no] WI 



0 10 20 30 40 50 60 70 60 90 


LATITUDE (® NORTH) 

I'igurc 1.3 Change in NC5 due to vertical transport by wave number 1. 


. . ■ 'T.-r?-; t - ^.s.-T :.rr :r 

N81-ft27F»4 

6 


2. ROCKET EXPERIMENTS 

1116 tasks of this group, under the direction of L. G. Smith, relate 
tc the definition of the scientific objectives of the Aeronomy Laboratory's 
rocket program; to the development and fabrication of instrumentation; to 
the scheduling and launch operations associated with individual and coor- 
dinated rocket experiments. This research is supported b) the National 
Aeronautics and Space Administration under Grant NGR 14-005-181. 

2.1 Introduation 

The Awronomy Laboratory participated in the Energy Budget Campaign v'ith 
throe payloads Loesignated Ell) on Taurus Orion rockets. The rockets were 
launched, as part of two salvoes, on 16 arid 30 November 1980, at ESRANGE 
near Kiruna, Sweden. The launch times and apogee altitudes are indicated 
in Table 2.1. Tlie throe payloads were identical except for the mass 
spectrometer, provided by the University of Bern, which measured positive 
ions in two flights and negative ions in the other, as indicated in Table 

2.1 by P and N, respectively. 

The principal objective of these flights is the study of the ionization 
and the ionizing sources (energetic particles) during periods of auroral 
activity. This is accomplished by measuring the concentration and temperature 
of the ionospheric (thermal) electrons, the ion composition, and the flux of 
energetic particles (discriminated by type of particle, spectrum and pitch- 
angle distribution). The payloads also included an experiment, provided by 
the University of Sussex, to study the high-speed fluctuations (bunching) 
in the energetic particle flux. We now give a general description of the 
payloads and the supporting systems followed by an account of the University 
of Illinois experiments and some preliminary results. 

The Energy Budget Campaign included other rocket launches, using the 
range at .\ndoya, Norway as well as that at ESRANGE; ground-based experiments; 
and some balloon-borne experiments. Also one pass of satellite S78-1 occurred 
within 3 minutes of the launch of Taurus Orion 33.009 on 30 November 1980. 

This will allow some interesting comparisons of energetic particle data 
from the rocket and the satellite. 

2.2 General Description 

The general arrangement of the Taurus Orion payloads is indicated in 
Figure 2.1. The diameter of the instrumentation section is 30 cm. At the 
front is an ejectable nose cone, protecting the input section of the mass 
spectrometer. Ejection takes place at an altitude of T 52 s (57 km, 
nominal altitude) and the cap of the mass spectrometer is ejected 2 s later 
(at 60 km) . 

Doors in the rear experiment section are also ejected at T + 52 s (57 km) 
and booms deployed at T + 60 s (67 km). The booms carry the electrodes of 
the probe experiment and the forward-looking energetic-particle detectors. 

The spin rate at boom deployment (and for the scientifically significant 


7 


Tabic 2 . 1 

Taurus Orion (lill payloads) 


ROCKIiT 

M.S.* 

LAUNCH CUT) 

SALVO 

■mnmi 

53.010 

P 

03:49:40 

16 November 1980 

B 

177.3 

33. on 

N 

04:15:00 

16 November 1980 

B 

179.2 

33.009 

P 

23:44:30 

30 November 1980 

A- 2 

171.7 


*Mass spectrometer: I’, positive ion; N, negative ion. 
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(BERN) 


TELEMETRY 
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(NASA) 
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(ILLINOIS) 


EXPERIMENTS 

(ILLINOIS) 


TELEMETRY 

(NASA) 

GYRO 

(NASA) 


RECOVERY SYSTEM 
(NASA) 


Figure 2.1 The configuration of the Taurus Orion (Ell) payloads 
showing the location of the instrumentation. 
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part of the flight) was specified to be 5.5 rps. The actual values (at apogee) 
for the three payloads range from 5.86 and 5.98 rps. The rockets showed very 
little preccssional motion. 

Hie payload includes a transition frustrum to adapt the 30-cm-diametor 
instrumentation section to the 36- cm- diameter recovery system. n»o second- | 

stage rocket ;Orion) is also 36-cm diameter. T1»e total length of the payload, | 

including the transition frustrum and the recovery system, Is 349 cm; the | 

total weight is 181 kg. Tlie nominal (predicted) apogee was 171 km, whid\ was I 

accomplished, and, for two flights, slightly exceeded (see Table 2.1). | 

'nie payload includes two independent S-band (about 2.3 GHz) FM/i'M tele- i 

metry systems. Tlie forward system supports the mass spectrometer experiment, i 

one tracking system (TRAPAT) and various monitors (bettery voltage, etc). j 

Hie rear system works with the Illinois/Sussex experiments, the gyroscope and 
a commutator which carries monitor signals (recovery system, etc). 

Hie three-axis gyroscope is located inside the transition frustrum. 'Il)e | 

yaw, pitch and roll signals are transmitted on channels 5, 6 and 7, respectively 
;uid the commutator on channel 8. ’llie output of a transverse magnetometer (to ' 

measure spin rate) is carried on channel 9, llie remaining channels (10 to 20, 
iuid 1) of the rear telemetry sj'stem carry signals from the experiments described 
below. 

I'wo tracking systems are used. Radar tracking by an MPS-36 is aided by 
a transponder located in the forward instrumentation section. ITte primary 
tracking system (TRAPAT) uses a PCM signal transmitted (at 500 MHz) to a 
receiver in the forward section of the payload. Pie information (3.9 kHz 
bit rate) is telemetered to the receiving station on channel 18 of the forward 
system. Riuige is obtained from the round-trip time delay. The azimuth and 
elevation angles are obtained from the S-band receiving antenna position. 

ihe recovery system severs the payload from the rocket motor at T ♦ 318 s 
(nominally 115 km, on descent). A second severance occurs at a nominal 
altitude of 6 km, releasing the main parachute. Tlie predicted impact occurs 
17.5 minutes after launch. 

Hie arriuigcmcnt of the rear experiment section containing the experiments 
of the Universities of Illinois :ind Sussex is shown in detail in ligure 2.2. 

Hie booms arc folded with the ejectable doors covering them. The booms are 
restrained by a cable which is cut after the doors have been ejected, n^e spin 
of the rocket causes the booms to deploy. At full extension the center-to- 
center distance of the probe electrodes is 179 cm. 

2. 3 MeasiO'enent of Thermal Electrons 

2.3.1 Introduction, llie probe experiment is used to measure the varia- f 

tion (fine structure) in electron concentration, to measure electron temperature * 
and in the measurement of vehicle potential. Ihc propagation experiment is 
used to obtain the electron concentration and the electron collision frequency. 

Botl) experiments arc well established in theory and in instrumentation. I 

I 
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Figure 2.2 The rear experiment section of the payloads of the 
Taurus Orion rockets. 
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2.3.2 DC/Langmuir' probe experiment, TTie electrodes osed in these pay- 

loads are 0.95 cm diameter and are coated with Aquadag (coliiodal graphite) 
to give a surface of uniform contact potential. The two electrodes are 
connected in parallel to the electronic circuits: this minimizes the spin 

modulation which appears in the probe signal. The voltage applied to the 
electrodes is fixed at 4V except that, at intervals of 2 s it is swept from 
-1.35V to 4V with a duration of 0.3 s. 

The probe current at fixed voltage is measured using a logarithmic 
electrometer and transmitted on channel 20. During the sweep period the 
rate-of-change of probe current is measured on the same electrometer. This 
allows the electron temperature to be obtained directly [Smith et al.^ 1968; 
Zirrxerman and Smithy 1980). This mode also allows vei'.icle potential to be 
obtained. 

ITie electrometer output is amplified during the period of fixed voltage, 
and transmitted on channel 18. The amplifier has a frequency range of 90 to 
2000 Hz [Klaua and Smith, 1978]. 'Hiis is the fine-structure experiment. 

2.3.3 Hadiouave propagation expeyd^ment. Ilic radiowavc propagation 
experiments used with rockets ai'e based on the techniques of dispersive 
Doppler and of differential absorption. The most successful techniques use 
the two characteristic modes of the wave whidi, in the ionosphere, have 
different refractive indices and absorption coefficients. In this context 
dispersive Doppler and laraday rotation are identical techniques. 

Tor the Lnergy Budget Campaign we chose to use a simpler technique than 
we had previously used, llus was done because of the availability of trans- 
mitters and antennas at LSRANGb. Ihe type of experiment has been described, 
for cxiunple, by Fvinae and Williame [1976] and Saoobscn atU Friedrich [1979]. 

Our pay load 4 each contain two independent propagation experiments: one 

operates at 3.8£3 MHz and the other at 7.835 MHz, llAe modulation outputs of 
the two receivers are transmitted on channels 14 and 15, and the AGC voltage 
on channels 10 and 11. 

'Hie dispersive Doppler and differential absorption techniques end when 
the extraordinary mode has been totally absorbed. In the three flights useful 
data are obtained between altitudes of 70 and 95 km. At greater altitudes 
some data may be recovered from the AGC signals since this can be related to 
absoiqition of the ordinary mode. 

2.3.4 Probe data ^v:ahJl^i8. A preliminary analysis has been done on 
the probe data from the three flights. Before a final analysis can be made, 
two computer programs used in reducing these data must be rewritten for a CDC 
Cyber 175. Tins is due to a change in the computer system used in digitizing 
the data at NASA (Wallops Flight Center). Also the original programs were 
written for an IBM machine. The University of Illinois now has a Cyber 175 
whidi can perform the same processing much faster (and cheaper) than the IBM. 

Shown in Figures 2.3 and 2.4 are the probe current profiles for the two 
Ir'onches during Salvo B: Taurus Orions 33.010 and 33.011. The altitude is 
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i0 !0 10'*^ 10'® 10'^ 10 ® 10 ® 10‘‘* 

OROBE CURRL'JT 

Figure 2.4 Probe current profile of 5.3.011 launched from liSRANGF 
16 November 1980 at 04:15:00 (UT). 
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in kilometers ami the current is in amps. Only the ascent portion of each 
flight is shown. The interrupted profile occurs because the sweeps (introduced 
to meas jre electron temperature and vehicle potential) have been removed. 

Very noticeable in Figure 2.4 is the sharp rise in current at 94 km, 
followed by a sharp decrease at 159 kra. It may also be noted that rocket 
33.010 (Figure 2.3) also crJiibits a sharp decrease in current at 139 km but 
no sharp rise at the lower altitude. 

The probe current profile for the launch in Salvo A-2 is shown in 
Figure 2.5. Again only the ascent portion is shown and the sweeps have been 
removed. Again there is a sharp increase in current followed by a sharp 
decrease. This time at 107 and 125 km, respectively. 

2.4 Encvgctio Particle Ex^H'riment 

2.4.1 IntvoductiaK. The detectors used in this experiment are silicon 
surface-barrier devices. Four detectors. Table 2.2, arc included in each 
payload: all have a sensitive area of 50 mm^, use a reverse bias of about 

lOOV and have a depletion la' er thickness of 300 urn. 

Two detectors are mounted on the booms (one on each boom; see Figure 
2.2) so that, with the booms deployed, the axis of the detector is 25 degrees 
from the payload spin axis. geometrical factor (GF) for the forward- 

looking detectors is 0.45 cm^ster; the large value was chosen to give a high 
count rate for the particle bunching experiment. The full width at half 
maximum (FWIIM) is 59 degrees. 

Tlie other two detectors are mounted inside the payload behind one of the 
door openings,^ iliese are oriented perpendicular to the spin axis. Tliey have 
GF = 0.057 cm"ster and FWllM » 18 degrees. Tlie narrower field of view gives 
better pitch-angle resolution. 

One of each pair of detectors has an increased thickness of aluminum on 
the front surface Tlie comparison of the flux seen by the pair of detectors 
allows some discrimination regarding the type of energetic particle (electron, 
proton, etc). Tins technique has been described by Fries ct at. [1979]. 

Pie general arrangement of the electronic circuits of this experiment 
is shown in Figure 2.6. The signal developed in the detector is a very small 
charge which is proportional to the energy of the particle (excluding the 
energy lost in traversing the metal layer). This signal is converted to a 
voltage in the low-noise pre-amplifier which has a capacitor as the feedback 
element. The pulse is now shaped (and noise is reduced) and further amplified 
so that a 128 keV particle gives a (negative) lOV-peak pulse. Pie pulse 
duration is about 3 us at this point. 

After amplification the pulses are processed in three separate systems. 

One is the particle-bundling experiment provided by M. P. Gough, University 
of Sussex. 'Pic other two systems are in parallel but are not redundant. The 
staircase generators indicate particles above a threshold energy at the instant 
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Figure 2.5 Probe current profile of 33.009 launched from LSRANGE 
30 November 1980 ar 23:44:30 (UTj, 
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Figure 2.6 ITie energetic particle experiment. 
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of occurrence. The microprocessor system, on the other hand, accumulates data 
over a period of about 5 seconds but sorts the data, before transmission, by 
energy and pitch angles. Thus the staircase output is useful for studios of 
pulsations whereas the microprocessor gives the energy rpectrum. Both give 
pitch-angle distributions. The telemetry channels used in this experiment are 
indicated in Figure 2.6. 

The accurate determination of particle flux for the side-looking detectors 
is difficult: the payload intercepts gyrating particles with pitch angles near 

90 degrees [VO 00 and Smith, 1977], 

2.1.2 Staircaoe 0 \j 0 tem. The pulses from each amplifier pass to a 
comparator set for -2V, giving a threshold energy of 25 keV. niis is a nominal 
value, the actual value is established by calibration with a radioactive source. 
Tlie pulses arc lengthened to 50 us in a onc-shot to establish a dead time and 
limit the maximum counting rate. An 8-bit counter (5 bits are used) and a 
digital-to-analog converter give the staircase output, one step for each pulse. 
After 32 counts, which is 5V output, the staircase returns to OV. tach of the 
four detectors has a staircase system. 

One advantage of this system is the extremely wide range of count rates 
that can be measured, 'Hie upper end of the range is limited on’/ by the in- 
formation bandwidth of the telemetry channel. Because of the onc-shot the 
maximum count rate that can be transmitted is 2 x 10^ s"^. Rates approaching 
this value must be corrected for the dead time. At this maximum rate the rate 
for the staircases themselves is less by a factor of 32 i.e. 625 s‘i. With the 
telemetry channels used there is a reduction in amplitude at this maximum rate 
but no loss of information on count rate. This reduction, and the strong 
pitch-angle dependence for the side-looking detectors is illustrated in Figure 
2.7. 


2.4.3 Hiavoprooeeaov Bjet&n. Tne Z80 microprocessor is preceded by a 
pulse-height analyzer. Four peak detectors. Figure 2.6, one for each particle 
detector, retain the pulse of largest magnitude. The four data channels are 
sampled sequentially at torvals of 50 ms. liach channel accumulates data for 
150 MS. It is then road and, in the remaining 50 ms of the cycle, is cleared. 
Ihis process is indicated in Figure 2.8(a). 

The exponents ul nature of typical particle energy spectrums make it 
appropriate to convert from the linear pulse-height measurement to a logarithmic 
one. This is simply accomplished using an erasable programmable read-only 
memory (EPROM) in which the inputs (linear data) are the address and the out- 
puts are the log data. 

The energy range is divided into 16 bins. The lowest contains only the 
counts from noise; the hightjst contains the counts from all energies greater 
than 128 keV. The remaining bins cover specific ranges of energy on a 
logarithmic scale. 

The energy spectrums obtained in this way are further subdivided by pitch 
angle. This is accomplished using the spin magnetometer as a reference. The 
spin period is divided into 32 equal intervals by a counter and the spectrums 
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Hisurc 2.7 The staircase output fron. u 

s ide- look IniJ detector and the 
magnetometer signal showing 
strong pitch -angle dependence. 
The photograph shows a partial 
overlap of a second scan. 
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accumulated by sector for a total time (^out 4.8 s) which is slightly longer 
than the time needed to transmit the data. 

Tlie data arc transmitted in serial-digital form. The energy spectrums 
for the four detectors are transmitted for each sector until all 32 sectors 
are completed. ITiis system is described in more detail in Bvaawell and 
:!rnth [1981]. 

llie flight data from the microprocessor requires special post-processing 
since it is transmitted as digital information over the FM/FM telemetry system. 
Tl\is is a form of frequency- shift keying (FSK). 

liacj byte of data was accompanied by stop and start bits so that a 
standard UART (Universal Asynchronous Receiver/Transmitter), such as the 
1013, ciin decode the serial format off the analog flight tape. 

As soon as the UART decodes a byte of data it delivers it as 8 parallel 
bits as its output port where an LSl/11 micro computer reads it in and stores 
the information on floppy disks. 

Hie information on the floppy disks is later transferred to a Cyber 175 
disk file for data iinalysis for FORTRAN programs. 

ilic high count rates encountered in the launches of the tnergy Budget 
Campaign, caused a problem for the pulse-height analyzer. As indicated in 
l-igurc 2.8, the smaller pulses are selectively ignored causing the apparent 
spectrum to be modified. If the pulses are assumed to be randomly occurring 
then a correction can be made to the data. 

Two additional features complicating the analysis of the particle data 
have been noticed. Tlie first, also a result of high coioiting rates, is that 
overflow occurs in the highest energy bin. T!\is is not serious but requires 
additional effort in the data analysis. The other, apparently unrelated to 
the high count rate, is that the azimuth reference appears to be displaced 
by about 60 degrees. Since this effect is also seen in the output of the 
staircase system (Figure 2.7) it is clearly not attributable to the micro- 
processor system. 

2.5 Future Plana 

The first post- flight meeting for the Energy Budget Campaign is to be 
held in Heidelberg, Germany in April 1981. 

Die analysis of data will continue to be a major effort in the next 
year. It is hoped to contribute papers to the COSPAR meeting of May/June 
1982 in Ottawa. Special sessions related to the Energy Budget Campaign are 
being organized by Dr. D. Offermann. 

Two rocket launches are planned for the summer of 1981. Tliese arc Nike 
Orions 31.014 and 31.015, instrumented to study energetic particle precipitation 
at middle latitudes and its effect on the nighttime ionosphere. 'Hie launeJaes 
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will take place near midnight from Wallops Island with the magnetic index 
K > 5. 

'nie experiments are similar to those used in the hnergy Budget Campaign 
(.see above) with the addition of a second probe experiment and a lunar 
aspect sensor. Both arc described in detail in Duvkin anJ Smith [1981], 

Tentative pljms for the more distant future include nighttime launches 
from Wliite Sjinds, late in 1982, re-using the linergy Budget Campaign payloads, 
n^e main scientific objective will now, however, be the yet unresolved 
question of the role of metal ions in maintaining the nighttime lower iono- 
sphere . 

In the following year (1983) it is noped to participate in a NASA 
launch operation at the magnetic equator using the existing site near Lima, 
Peru. 
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3. PARTIAL RLFLECTION 

This group is involved in the study of the d region of the ionosphere 
by the techni(|ue of partial reflection, lliis includes regular measurements 
of electron concentration by the differential-absorption technique for use 
in studies of region variability, further work on the differential-phase 
technique, improving measurement techniques for the lower V region, and 
part ial - reflect ion drifts measurements. Tlie research is supported by the 
National Aeronautics and Space Administration under grant NCR 14-005-181. 

3. 1 I'rbana Meat^urerjenta 

K. M. Wei land has collected differential-absorption electron concentration 
and drifts wind data on a daily basis near local noon at the Aeronomy 
Laboratory Field Station near Urbana (geographic coordinates: 40® 10' 

10" N, 88“ 09' 36" W) , from December 1980 through mid- February 1981. Data are 
digitized at 1.5-km intervals in the altitude range of 60 to 90 km. Files 
of 512 samples (or 3.4 minutes) are collected with 0 dB, 10 dB, and 25 dB 
of attenuation ahead of the receiver. Tlie lowest attenuation file which doei 
not saturate the receiver at a given altitude is retained for further processing, 
live of these groups of unsaturated composite files (for a total of 15 files) 
are collected. Tlie collection time is about one hour. The median of the 
five unsaturatod values of electron concentration at each altitude is found to 
obtain a final e lectron-coneent rat ion profile. Ilic collection routine has 
been further described in ;ut earlier report utiJ Houhi 1 1 ^ 1973]. 

During the past winter the data collection time was doubled. A partial- 
reflection drifts wind measurement was first made, followed by an hour of 
electron-concentration data collection, a winds measurement, a second hour 
of electron-concentration data, and ending with a winds measurement. Winds 
were measured at 3-km intervals from 75 to 90 km. The collection procedure 
was modif.ed to fill one half hour during the approximately 25 days that the 
Urbana meteor radar was operational during this winter. 

Figures 3.1 through 3.3 show the electron-concentration data obtained 
during the months of December 1980 through February 1981. Data plotted 
includes the electron zoncent rat ion at 72, 76.5 and 81 km, and the 
ratio at 81 km. 

The daily winds data are being compared with the daily electron-concen- 
tration data to study the meteorological causes of the winter anomaly. By 
avo'-aging the .hree wind profiles per day it was hoped that the effects of 
short-period fluctuations in the wind would be removed to provide a wind 
estimate more closely approximating the prevailing wind for that day. 

I’reliminary comparisons of these winds data with electron-concentration data 
liave not shown very correlation, indicating that the noon winds are not 

a good predictor of the 24-hour prevailing winds. 


The Urbana meteor radar system was operating for about 25 days during this 
l^ast winter. The meridional prevailing wind is being compared with the 
electron concentration measured on those days. Figure 3.4 is a plot of the 
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meridional wind component and 81 km (inversely proportional to the total 

electron content below 81 km). It imist be noted that in this preliminary 
figure the winds data plotted are from a 24-hour period centered at noon on 
the day shown. For this conqiarison the winds for a 24-hour period centered 
at midnight the day bo<^ore arc needed, so when viewing Figure 3.4 the winds 
should bo shifted to the right by one day. The correlation between winds 
and electron concentration shown in Figure 3.4 is quite good. 

3.2 Future Syetem Development 

Investigation is being made into emulating the partial-reflection drifts 
experiment on an Apple microcomputer. Preliminary work has been done by 
R. Ruggerio into the various methods which can be used to implemen: this 
system. 

The present method is that of the spaci)d antenna. This method involves 
the analysis of the amplitude information received from four spaced antennas. 
By computing the auto- and cross-correlation of the returned signals, the 
wind velocity is computed using the full correlation analysis. 

The other method involves using coherent integration. This requires 
the addition of phase detectors to the present system in order to extract 
the in-phase and quadrature components so that a complex correlation is formed 
for the determination of the drift velocity. 

The implementation of the microbastd system will allow simultaneous 
measurements with meteor radar or coherent scatter which wili enable accurate 
correlative measurements. 

The analog-to-digi tal interface for this system is under construction. 
This interface will contain two channels making it compatible for either the 
amplitude or phase method. Height resolution of 1.5 km is expected from the 
interface which is identical to the present system. 

If the amplitude method is used the second channel in the interface can 
be used as a backup circuit or to implement dual processing by using two 
microcomputers simultaneously. Dual processing may be implemented if the 
number of altitudes desired is more than a single microcomputer can handle. 
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4. METOOR RADAR 

Experimental and theoretical studies of the dynamics of the meteor region 
(75-105 kmj are being conducted using a meteor- radar system. Ibis research 
is supported in part by the National Science Foundation under grant ATM 79- 
09982; in part by the National Oceanic and Atmospheric Administration under 
grant NA-80-RAD-00016; and in part by the National Aeronautics and Space 
Administration under grant NCR 14-005-181. 

4.1 Meteor Radar Obeervatione 

During the periu.; October 1, 1980 to March 31, 1981, the radar has been 
operated approximately 700 hours, covering the periods: 

November 4, 1980 
December 2-5, 1980 
December 16-19, 1980 
January 20-23, 1981 
i ebruary 3-8, 1981 
February 17-20, 1981 

The data have been fully analyzed to yield the north/south and east/west 
horizontal wind components, and curve-fitting routines are used to yield the 
daily mean winds, and the diurnal and semidiurnal tides. 

4.2 LuLiV Correlative Ctudy 

As discussed in Aeronoray Report No. 80-2 [Eduarda, 1980], near-simultaneous 
meas.irements of the wind field and sodium density were obtained on January 19, 
1980, March 15 and March 22, 1980 in order to study the effect of dynamics on 
the structure of the sodium layer. The results have been further analyzed 
and are summarized as follows: Ibe east/west wind field has been reconstructed 

using rhe mean ♦ diurnal ♦ semidiurnal components, and compared with the sodium 
density profile at the same time. If the wind-shear mechanism for the down- 
ward sweeping of ionization is dominant, then one should see peaks in sodium 
density at roughly the same altitude where there is a wind mode, with west- 
ward winds above eastward winds. 

On March 15, there is a node at 90 km with thj requisite westward above 
eastward winds (Figure 4.1a). At the same time, the principal peak in the 
sodium density is also at about 90 km (Figure 4.1b). The profile shows very 
little wave activity. Similarly, the vertical wind profile on March 22 shows 
a node at about 94 km with westward above eastward (Figure 4.1c). The 
principal peak in the Na profile is at approximately 93 km, again showing 
good correlation with the tides (Figure 4. Id). More wave activity is present 
in the sodium profile, and a spectral decomposition of the north-south wind 
shows that the high-frequency components are strong relative to the diurnal 
and semidiurnal tide that day (Figure 4.2), as compared to the data of March 
14-15 (Figure 4.3). If this correlation, between little wave activity in the 
Na profile and little high-frequency wave activity in the horizontal winds, 
carries sufficient weight, then we are in a position to separately identify 


altitude Uffl) 
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MANCH »t, mo. 0049 CST 



Figure 4.1 (a) Mean ♦ diurnal + semidiurnal winds at nir^O hrs CST on March 15, 

1980; (b) sodium density profile; fc) mean ♦ diurnal semidiurnal 
winds at 0049 hrs CST on March 22, 1980; (d) sodium density profile. 












spectruni for March 14-15, 1980. 
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tidal and gravity wave effects, since the wind shear mechiinism does not apply 
to north-south ion drifts. Tlie high-frequency waves in the periodograms are 
definitely gravity waves. 

l-inally, the data from January 18-19 show a node in the wind profile 
at roughly 93 km (Figure 4.4a). At the same time, the sodium density profile 
shows a peak at 90 km (Figure 4.4b). Tliis discrepancy is unresolved at the 
present time. However, the spatial periodogram reveals a vortical wavelength 
of about 9.1 km, whereas the analysis of the vertical structure of the semi- 
diunial tide shows a wavelength of about 10.3 km, with downward phase propaga- 
tion. Tliis motion is evident in the lidar data, and it is plausible to suggest 
that the Na layer is indeed tracing tidal motion via the wind shear mechanism, 
in addition to tracing gravity wave density perturbations. Periodognims are 
noisy estimators of spectral information, and it is not unreasonable to suggest 
that the 9.1 and 10.3 km wavelengths correspond to the same phenomenon. 

Hic results of these correlative studies, albeit modest, do seem to 
point in the right direction, and further analysis of meteor radar winds may 
determine which features of the sodium profile are due to gravity waves or 
due to tides. 

4.3 Ir:pvoih^''iento 

An improved data collection system is presently being designed for the 
Urbana meteor radar experiment. The present system uses a Digital Fquipment 
Corporation PDP-15 for real-time data collection and processing. nic new 
system features two Apple 11 microcomputers, one for raw data collection 
and the other for preliminary processing. 

Hie new system will provide the opportunity for simultaneous data 
collection with other experiments and will be more reliable than the old 
system. At the present time, correlative studies between exjieriments must 
rel)' upon interspersed data collection due to equipment conflicts. Die new 
setup allows for simultaneous meteor-radar and partial -ref lection data 
collection along with other experiments in the future. Die reliability of 
the system will improve because of its modular design. Any component that 
fails can be quickly replaced. 

Most of the work to date has gone into the investigation of the various 
languages available for the Apple II. Tlie FORTH language provides low-level 
control of mcmoi-y and is compatablc with machine language. It seems to be 
the best choice for the raw data collection. BASIC and i’ASCAL are higher 
level Ian .uages vvhich arc easier to read, understand and maintain. Tlie FORTH 
language is a faster processing language and it has now been extended to 
compute floating point operations. The processing language has not been 
selected . 

Some of the hardware choices for the new system have been explored. A 
sample-and-hold array has been fabricated and bench tested, but final selection 
of the interfacing hardware must wait until the collection software has been 
written and is available for testing. 


ALTITUDE (km) ALTITUDE (km) 








4.4 The Pokev Flat, Alaeka MST Radar as a Meteor Radar 


Uuring the past eight months, an echo selection criterion has been 
developed to distinguish between cdioes due to meteors from those due to 
turbulent mixing (regular eclioes) that are received on the Poker Flat, 

Alaska MST radar. Several interesting results have been found. The dai ly 
variation of edio rates meeting this criterion can be partially explained 
by the geometry of the earth's orbital path intersecting the meteor trajec- 
tories and the geometry of the two antenna directions. The height profile 
of meteors agrees well with what is expected from the ablation process of 
meteors entering the atmosphere. Winds calculated from the doppler shift 
of the trails are not completely in agreement with the winds determined from 
the regular echo data. More work needs to be done in this area. 

The selection criterion is schematically shown in Figure 4.5. A meteor 
edio must not appear on both receivers at (tp* ho), it also must not appear 
at (to. h()3 and (tQ, ho^U and (tg, hQ-l). Nor can it appear at (to, ho) 

.md (tQ-1. lip) • Pinally, it cannot appear at more than two 

of the eight remaining positions shown in the figure. Since each grid 
point represents a four-minute average of data, this criterion is probably 
too exclusive; more work will be done on relaxing the assumptions in the 
cri terion. 

Based on this criterion, histograjns of daily meteor rates and height 
distributions were determined, Tlie height distribution of the meteor echoes 
is shown in Figure 4.{). Tliis profile can be explained in terms of the 
physical characteristics of meteors. As a meteor enters the atmosphere, 
heating and ablation begins. Depending on the mass and velocity of the 
meteors, these processes begin at 100-120 km where the air density is large 
enough to generate friction. As the meteor penetrates further, the expo- 
nentially increasing air density leads to an increase in ablation and 
ionization imtil the ionization reaches a maximum at roughly 92 km. At 
lower altitudes, the ionization falls off rapidly as the meteor diminishes 
in size. Die small thickness of the meteor region (80-110 kmj is the result 
of the rapid change in air density. 

llic daily meteor count for each receiver is shown in Figure 4.7. 'Hie 
hourly distribution is partially consistent with the earth's geometry and 
with the directions of the two antennas. The daily variation is not con- 
sistent with what is observed with the Urbana meteor radar system whidi 
peaks at 0500 and 1100 CST. Tliis may be due to errors in the selection 
criterion and/or poor statistics due to the low echo rate on the Poker Flat 
MST radar. Ihe Urbana meteor radar has a maximum echo rate of approximately 
500 !ir~^ compared to 4 hr“^ from the Poker Flat radar. IVhen the MST radar 
is in full operation, the edio rate should improve by a factor of at least 
100. liven with the current low echo rate, by looking at several days, a 
measurement of the prevailing wind and diurnal and semidiurnal components 
is possible. 


As meteors enter the upper atmosphere, they vaporize forming trails of 
ionized particles. Tlie trails diffuse rapidly but during their brief 
existence they reflect radio signals. Because the neutral-ion collision 
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frequency is so high (200 collisions s"^) the ionized trail drifts with 
the neutral wind. Consequently, the Doppler shift of the return signal 
gives a measure of the wind motion. 

'Hie mean wind for days 160-196, 1979 and days 173-205, 1980 is shown in 
Figure 4.8. Also shown is the mean wind determined from the regular meso- 
spheric echoes at Poker Flat and the mean meridional wind determined from 
the Urbana, Illinois meteor radar system during the same period. First, note 
that there is intcrannual variability in the mean wind. Tlie eastward wind 
determined from the meteor echoes has the same trend as the wind determined 
from the regular eclioes, however, it is approximately 20 ms"^ too low in 
comparison. This is true for both the 1979 and 1980 data and needs to be 
investigated further. Hie northward wind from the meteor echoes is also low 
when compared with the regular echo wind but the differences are not as great 
as those found in the eastward wind. A comparison with the Urbaiia northward 
wind shows that the mean wind is weaker at lower latitudes. Hus is in agree- 
ment with theoretical studies of the meait circulation in the middle atmosphere 
as discussed in Schoebcvl and dtvobel [1978]. 

Hie north and cast diurnal tides for the same period are shown in Figure 
4.9. nie eastward diunial tide from the meteor echo data has an amplitude 
that is low compared to that calculated from regular echo data. Also, the 
phase structure is scattered. Hiis may be due to the daily variability in 
the 24 hour period wave. Tlic northward diurnal tide compares more favorably 
both in amplitude and phase although values are still somewhat lower. The 
Urbana diurnal tide showed very little daily variability during this time 
interval based on daily Fourier jinalyses. Due to the low meteor echo rate 
at this time, daily Fourier analyses on the Poker Flat data cannot be cal- 
culated to compare with the Urbana results. 

Hie eastward and northward semidiurnal tides are shown in Figure 4.10. 
Again, the eastward amplitudes are weaker from the meteor data but the phase 
structure compares fairly well. The eastward phase structure indicates an 
upward propagating energy wave below 88 km and downward energy propagation 
above 88 km. Tlie regular echo data indicates only an upward propagation but 
it must be remembered that the regular echo data is approaching its upper 
height limit at that point. Tlie northward tides compare more favorably. 

The differences between the Urbana and Alaska tides are seen in the phase 
structures here. The Urbana semidiurnal tide has a much longer vertical 
wavelength than the Alaska tide since its phase profile is nearly vertical. 

Tlie observed tidal structure has been compared with recent theoretical 
calculations from Forbes [1981]. The modeled eastward diurnal tide shows a 
strong upward phase propagation. ITie Poker Flat regular echo data has this 
same trend, but not as strong a phase tilt. Both the regular echo and 
meteor echo northward diurnal phase agree well with the modeled results 
indicating a constant phase height profile near 1200 hrs. The Urbana north- 
ward diurnal tide also has constant phase with increasing he* get which dots 
not agree with Forbes’ calculations which show a strong dov.nward phase 
propagation. 


EASTWARD NORTHWARD 



-20 -10 0 10 20 -20 -10 0 


AMPLITUDE (mf') 

POKER FLAT, ZnmPOKER FLAT, REGULAR ECHOES URBANA 

METEORS 

= 1979, DAYS 160-196 = = =1980, DAYS 173-205 


4.8 Mean eastward and northward winds. 
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SEMIDIURNAL TIDE 
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Fi gure 4.10b 


Northward semidiurnal tide. 
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The semidiurnal eastward tide agrees well with the calculations at the 
lower levels, however, the modeled tide does not reverse phase direction 
above 88 km. The phase structure of the Poker Flat northward semidiurnal 
tide also is in good agreement with the theoretical calculations. The Urbana 
northward tides do not have the downward phase propagation as indicated in 
Forbes' calculations. 
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5. COHERENT- AND INCOHERENT- SCATTER RADAR 

Experimental and theoretical studies of the ionized regions of the 
atmosphere are being conducted using high-power backscatter radars at VHP 
and UHF frequencies. The first year (1978) of data from the Urbana coherent- 
scatter radar have been analyzed by 0. R^yrvik, K. Gibbs and S. A. Bowhill. 

This research is supported by the National Science Foundation under Grant 
ATM 78-21765 and ATM 78-15224 and by the National Aeronautics and Space 
Administration under grant NSG 7506. 

5 . 1 /irrai/ Prooeeeor Reeearah 

5.1.1 Introduction. The needs of computer users are demanding faster 
and more powerful computers. Faster computers may be obtained by either 
improving the raw circuit speed or by using the same circuits in a more 
efficient architecture. Unlimited improvements in circuit speed are not 
possible due to fundamental physical constants such as the speed of light. 
Therefore new approaches to computer organizations must be found if projected 
demands of computer users are to be met. 

Only one particular organization is dealt with in this study by R. 0. 
Colbert: the array processor. In this form of processor, high throughput is 
achieved by introducing parallelism, that is to say, several processors per- 
forming nearly identical operations. In this work, the array architecture is 
examined on a new approach to the design of an array processor is proposed in 
order to take advantage of the recent advent of low-cost, high-performance 
microprocessors . 

I Iliac IV is taken here to be the conventional array processor (Figure 
5.1). In an array processor, only the arithmetic units and memories are 
replicated and one single control unit (CU) drives the array of arithmetic 
units. For certain arithmetic units, the operation may have to be modified 
or suspended based on the contents of a gate or mark register in each processor. 
Hence, a minimum amount of control is kept local and each arithmetic unit plus 
its local control is called a processing element (PE). 

Greater flexibility is obtained through the use of local indexing, mode 
control and routing. Local indexing allows each address broadcast by the 
CU to be modified by each PE. Mode control allows each instruction to be 
modified by each PE. Routing is defined as the method by which PEj may obtain 
an operand which is stored in PEMj (i j). 

5.1.2 Issues of the study. Several issues related to the array 
architecture, such as processor* memory interconnection, interprocessor 
communications, software requirements, system expandability and fault 
tolerance of system. This study examines the results when the processor- 
memory interconnection consists only of a single address bus and a single 
data bus. This research also explores the problem V'fhich results when the 
only form of communication between two PEs is through the stored memory 
(SM). 'Phis work considers the possibility of recoding uniprocessor programs 
into code which contains nc context-dependent branches. That is, all PEs 
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must branch the same way when they encounter a branch instruction else the 
branch must be removed and that segment of code recorded. This work examines 
the expandability of on array processor which executes only Context- Independent 
Code and the possible fault tolerance and fault detection capabilities of 
such a system. 

5.2 Approach to a New Array Proaeaaor 

Context-dependent branches reduce on array system’s throughput 
significantly. With no context-dependent branches (assuming little or no 
memory contention), the throughput of n processors is n times that of a 
single processor system. A progriya containing context-dependent branches 
requires that some of the PEs must wait while the other execute the branch 
of code they desire and then the groups reverse roles forcing the other PEs 
to wait while the first PEs execute their branch of code. The worst possible 
case for an array processor is the case of nested-context-dependent branching. 

5.2.1 Context independent code and its implementation. A new concept 
in the generation of array processors is introduced here. This is the 
concept of Context Independent Code (CIC). Tlie principle behind CiC is the 
elemination of context-dependent branches. One notes that CIC eleminates 
all context-dependent branches not all condi tioned branches from the program. 
Context-dependent branches are those conditioned branches for which the data 
or condition may be different for different PEs. Conditioned branches which 
cause the program to loop a certain number of times are not context-dependent 
as all PEs branch the sajne way every time. Context-dependent branches must 
be recoded so that the program appears not to branch at all. For instance, 
the usual algorithm for multiplication shifts the multiplicand and tests eacli 
bit of the multiplier. If the multiplier bit is 1, the multiplicand is 
added to the partial product; if the bit is 0, the program skips the add 
instruction. For the conventional array processor, all PEs whose bit is 1 
execute the odd instruction while the rest of the PEs are turned off. 

This program written in CIC causes all of the PEs to execute the odd 
instruction, the difference being that those PEs whose bit is 1 add the 
multiplicand to the partial product while those PEs whose bit is 0 add zero 
to the partial proi'uct. 

One might argue that adding zero is no better than turning a PE off. 
However, ti e CJC approach accomplishes two things which the previous approaches 
have not done. 

1) All PEs are constantly synchronized in lock step mode. 

2) The software is tailored to suit the array processor rather than 
tailoring the processor to suit sequential software. 

The single restraint of CIC is that all PEs must execute the same 
instruction at the same time. This means that the operands of the instruction 
determine which brands the PE is actually executing. Several examples of CIC 
are included in the study, among these are 8*8 multiplication, 8-bit 
Magnitude of Twos -Complement Numbers, 16/8-bit Binary Division, 32-bit 
Accumulation and 32 x 32 bit Multiplication. 
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5.3 Hardijave Aepccfta of tho Super^-Sb Multi-Miavopr*ooeeeor Syatm 

After careful analysis of the available, single-chip, 8-bit microprocessors, 
the MOS Technology 6502 microprocessor was selected for the array processor. 

Selecting the 6502 allowed the use of a versatile microcomputer known as 
the APPLli 11. The Al’PLt II provides users with eight peripheral connectors 
along the back of the APPLP's main board. Iliese slots allow the user to 
plug in proto-boards containing user-designed circuits. The APPLE II has 
available a supervising program called the APPLE System Monitor. Use of the 
APPLE II relieved the burden of writing a monitor for the array processor. 

Ihe APPLE II also provides the user with a mini-assembler which allows one 
to type programs into the APPLE in assembly level language. For these and 
other reasons, the APPLE II is an excellent choice for implementing an array 
processor. 

5.3.1 Avoh: t rotuih' of the ovcx'all ayctrr". The Super-65 consists of 
i'our 6502 microprocessors, one being the APPLE II 6502, with three others. 

Each 6502 is given a private RAM of dimension IK by 8 bits. 1'he rest of the 
6502' s memory (the upper 65K) resides in the APPLE II. Iliis is both RAM 
and ROM. llic combined RAM ami ROM arc jointly designated as Stored Memory 
(.SM). Hk' ardiitecture of the Super-65 (Figure 5.2) allows eacli processor 
to access its private me.iiory at any time but only one processor may access 
the Stored Memory Address Bus at a time. All processors access the SM 
Data Bus on a Read. However, for obvious reasons, only one processor may 
access the SM Data Bus during a Write operation. The processor having 
control of the Stored Address Bus and the Stored Data Bus during a Write 
operation is called the Controlling Processor (CP). The CP has control 
over all stored resources. A convenient feature of the architecture is 
that it requires no modification of the APPLE II. Tlie reason is that the 
APPLE 11 6502 is removed from its socket and placed on a peripheral board. 
However, a 40-pin conductor connecting the 6502 to the empty socket allows 
the Al’PLE to operate as if the 6502 were actually still in the original 
socket . 

Hie APPLE II 6502 is in control when the APPLE II is powered up, just as 
it would normally be. In order to designate another of the processors as 
the Cl’, one simply addresses the memory range $Cnxx, where n equal the 
processor number (1, 5, 5, or 7). Seven is the number given to the APPLE II 
6502. IVlten one accesses Cnxx i,x « don't care), the APPLE II automatically 
activates the I/O Select line on the n th peripheral slot. ITiis signal 
is used to give the processor in slot n control of the stored buses. When 
one accesses $C4>n4> (n = processor number 8), the APPLE II automatically 
activates the Device Select line on the n th peripheral slot. Tins signal 
is used to stop the clock of the n th processor. Wlien one accesses $CFFF, 
the APPLE II automatically activates the I/O Select line on all of the 
peripheral slots. Tlus signal is used to resynchronize and restart all of 
the processors simultaneously. 

As the execution of CIC programs requires that all of the PEs have the 
same address at the same time, a minimal amount of fault detection was 
implemented by simply comparing the address on eacli of the private address 



Figure 5.2 Super-65 system block diagram. 
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buses to the address on the stored Address Bus and lighting on LED if they 
did not match. 


Finally provisions has been made for the inclusion of a 6522 Versatile 
Interface Adapter (see Figure 5.3) which contains two 8-bit parallel ports, 
on each of the processor boards. 

5.3.2 Sumriary. The main objective of this work was to assort that 
recoding of microprocsjssor programs into CIC programs is feasible for an 
important set of applications. Iliis objective was achieved by implementing 
a four-processor array and recoding several programs for it. The programs 
were divided into two categories, independent and dependent data handling, 
llie first catagoj 7 allowed ead\ processor to work on a separate set of 
data so that no processor required any results from any other processor 
in order to complete its task. The second category required each of the 
processors to work on a subset of the entire problem. Tliis merjit that the 
processors needed to conununicate at various times in order to complete their 
task. It was found that for the set of independent-data recoded programs 
the throughput of an array of n processor would range from 0.45 n to n times 
the throughput a single processor. For the set of dependent data recoded 
programs the n processjr array would achieve a throughput of approximately 
0.3 n to 0.6 n (see Tal le 5.1) as long as the carry propagation delay is 
minimal. For the architecture chosen, their delay will not be minimal 
as n gets very large and this is one weakness of this architecture. One can 
sec then that so long as carry propagation delay iind the time required to 
store final results in Stored Memory are insignificant, one may achieve 
roughly a linear increase in throughput by adding more processors. 'Hiis is, 
of course, after the initial penalty for the recorded program has been taken 
into account. In theory, then, one may achieve arrays of thousands of 
processors each of which adds on incremental throughput to the overall system. 

Hiere are several areas for further researclu Among them is the 
pursuit of the ideal microprocessor solution of the carry propagation delay 
problem and of the problem of storing to Stored Memory and the recoding of 
more sample CIC programs to obtain a better understanding of the potential 
for and against CIC programs. 


5.4 


v.v 


'pjet'n ufid APPLh' 


Hie final stage of post-processing for the Urbana coherent-scatter radar 
is presently done on a Hewlett-Packard 9830A desktop computer with pen-plotter 
iuid paper tape reader. Brief descriptions of the post processing are contained 
in J:}’bs .u;.;’ ' ' [1979] and FAxu'ds [1979]. 


K. Gibbs has built an interface between the HP9830A computer and an 
APPLE II plus microcomputer. There are several reasons for having a 
bidirectional link between these two machines: 1) the pen-plotter attached 

to the 1IP9850A can be used from the APPLE II computer; 2) the APPLE II, 
which is a much faster madiine, can be used for more complex post-processing; 
3) the data previously collected can be transferred to APPLE 11 disk storage 
to be compatible with an APPLE II based data-col lection system under 
development. The design of the interface is discussed on the next page. 
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Table 5.1 Recoding factors for the sample program. 



Independent Data 

Dependent Data 

3-Bit Magnitude 

0.33 


8 X 8-Bit Multiply 

0.53 


16/8-Bil Divide 

0.45 


32-Bit Accumulator 

1 .00 

0.60 

32 32-Bit Multiply 

0.59 

0.30 
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The HP9830A section of the interface ia a standard HP11202A 1/0 TTL 
interface card. The card is bidirectional with separate lines for input and 
output. The inputs are terminated with an active pulliqp and the output lines 
are driven open -col lector. All signals on the HP lines are inverted. A 
handshake control pair is also provided on the card. 

The APPLE II section of the interface is shown in Figure 5.4. It con- 
sists of a 6522 VIA chip and several TTL chips tc adapt the 6522 to the HP 
line conditions and the APPLE 11 bus timing. Port B of the 6522 drives open- 
collector inverters to transfer to the HP9830A. Lines from the HP9830A are 
activel/ terminated, driving port A of the 6522 through inverters. The 
74LS74 idjusts the timing signals from the APPLE 11 to properly run the 6522. 
Finally, the control pins are wired for a handshake mode for either 
direction of transfer. The interrupt flags for both parts are raised to- 
gether however so that switching data directions requires clearing the flags 
before starting. 

5.5 Conetmation of VHP Radar Preprooeeaor 

During this period much of the preprocessor was constructed and debugged, 
llie hardware consists of four 4V * 9" vector boards with standard wire 
wrapping techniques used for interconnections on each board. The four boards 
are contained in a standard equipment rack chassis allowing interboard 
communication via a wired backplane of the four boards, Board *1 contains 
the preprocessor memory (4K RAM), the memory address register (MAR) and the 
arithmetric logic unit (ALU). Board #2 holds the control store (two 
32 X 8 bit ROMs), instruction decoders and control logic. Additional 
control logic resides on Board #3. Board #4 contains the sample and holds, 
analog- to-digital converters and a pair of latches. 

Two valuable debugging and troubleshooting aids were included in the 
preprocessor design, single step capability and external access to the 
instruction bus. The single step circuitry disables the system clock and 
allows the user to issue individual clock pulses. The preprocessor executes 
one instruction for each clock cycle hence by single stepping the pre- 
processor can be examined before and aftrr an instruction execution and the 
change in logic state determined. 

External access to the instruction bus enables the user to issue any 
desired instruction to the preprocessor. Combined with single stepping the 
user is able to test any sequence of instructions desired observing the 
change in logic state for each instruction. This capability was used to test 
the microcode which will control the preprocessor operation. At this time 
the preprocessor correctly single steps through the entire microcode program. 
Tl\e next step is to program the control store PROMs and run the preprocessor 
from its 10 MHz system clock. This work was done by F. Zendt. 

5.6 Neu Coherent-Soatter Antenna 

The specifications and preliminary design configurations for a new 
coherent-scatter phased array antenna have been reported in Eduard^ [1980] . 
The proposed antenna shall have an aperture of 40.92 MHz over a 1 MHz band- 
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width. The array will bo steerable in two orthogonal planes to discrete 
angular directions of 30 degrees from the vertical, in addition to a 
vortical pointing direction. 

A proposed realization is to space half-wave dipoles on a square grid 
of one half wavelength spacing, a quarter wavelength above ground plane. 

To scon the beam in one plane, progressively phased excitation currents of 
0. 90, 180, and 270 degrees are applied along each row of elements Like- 
wise. applying the progressively phased excitation along the columns scans 
the beam in the orthogonal plane. Feeding all elements in phase produces 
a vertical beam. This arrangement results in every fourth element, in each 
row or column, being feed in identical phase for each pointing direction, 
llio inphaso elements will be hardwired together into 16 interlaced arrays. 
ITierefore, emanating from the antenna are 16 ports which scan the bean when 
excited with the appropriate phases. 

To evaluate the performance of such a system, a scale model has been 
built for antenna pattern studies by D. Tanner and J. Loane. To cover the 
required aperture the full scale realization will contain in excess of seven 
thousand dipoles. The model, however, has been limited to a small section of 
64 dipoles (8 x 8) due to practicality and constraints of the available 
antenna range and equipment. The scaled wavelength of the model is exactly 
6 inches corresponding to a frequency of approximately 1.97 Qlz. TTie unsealed 
wavelength is 24 feet; thus the frequency scaling factor is about 48. 

The model is fed through 50 Q coax cables of matched electrical length. 
Thirty-two half-wave baluns were etched on Rexolite PC board to convert the 
unbalanced feeders to a balanced 200 Q source. A pair of dipoles of 2 
wavelength separation is fed in parallel from each balun through equal trans- 
forming lengths (i.e. and odd number of quarter wavelengths) of balanced 
transmission line, llie proper phasing of the input cables is achieved with 
a network of power dividers and quadriturc and 180 degree hybrids. 

5.7 Midlatitude Mcaoephevia Scattered Power and Winds 

VIIF radars have been used successfully to observe the dynamics of the 
mesosphere at the equator [Woodman and Guillen, 1974; Raetogi and Bowhill, 
1975a, b; Cou/itryman and Bowhill, 1979], at midlatitudes [Miller et at., 1978; 
R^ttger et at., 1979] and at high latitudes [Balaley and Gage, 1980]. Al- 
though VHF radars can detect incoherently scattered signals from thermal 
fluctuations in the electron density, most VHF radar observations from below 
100 km rely on coherently-scattered signals from irregularities in the 
refractive index caused by turbulent advection of ionization gradients 
provided by the continuous increase with altitude of the mixing ratio of 
electrons between 60 and 90 km. Coherent-scatter radars are able to obtain 
1 ine-of-sight velocities at these altitudes, and thus are important tools 
in the study of dynamics of the mesosphere. 

A coherent -scatter radar has been operated by the University of Illinois 
at Urbana since the spring of 1978 [Miller et at., 1978; Gibbe and Bowhill, 
1979], and data have been collected during daylight hours on a number of days. 
This section by 0. R^yrvik describes the spatial and temporal variations in the 
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scattered power, and their relation to atmospheric motions of different 
amplitude and period in the mesosphere. 

5.7,1 Piivev and aot't'ePita’^^n *ime. VHP scattered power returns from 
the mesosphere over Urbana vary substantially as a function of both altitude 
and time, llie most prominent feature of the power profile is the consistent 
occurrence of peaks, or layers, of high scattering cross section. ITie 
scattered power in these layers are typically 10 dB above the reference level, 
but ciui reach as much as 25 db above that level (Figure 5.5). Die sporadic 
enhancements in scattered power must represent a local increase In the 3.7- 
meter-wavelength component of the turbulent spectrum. The tendency of these 
"whitecaps" to occur in well-defined layers is illustrated by Figure 5.6, 
which shows hourly 50** and 90"* power profiles for four hours. These layers 
of high scattering cross section may remain stationary in altitude for 
several hours, as is seen in the peak at 80 km altitude. Ihere may be small 
but significant vertical shifts in the scattering layer within one hour, as 
can be seen in Figure 5.5. These apparent vortical motions are generally 
small, one or two kilometers or less, and the shift may take pla:e within a 
minute. 

It is necessary to be cautious in interpreting these apparent vertical 
motions of the layer. Ihey certainly cannot be vertical mass motion of the 
air itself; they may however, correspond to turbulence associated with either 

(a) a vertically migrating shear; or 

(bj a sloping shear layer which is being carried past the vertical 
radar beam by a horizontal wind. 

Since wc cannot easily resolve between these two possibilities, we will refer 
to the phenomenon as "vertical motion" even though no actual vertical motion 
may be involved. 

Since we know that tidal oscillations and long-period gravity waves 
produce downward-propagating velocity pi'ofiles, it is tempting to explain 
the vertical motions of the layers on this basis. However, seldom is such 
a continuous downward propagation seen. Table 5.2 shows the number of days 
on which different types of behavior of the scattered-power profile were 
observed. It is seen that continuously descending layers were found on only 
1()1) of the fall 1978 data and none at all in the spring data set. 

The best example of a continuously descending scattering layer is seen 
in Figure 5.7. A strong stable layer was present at the start of the 
observation at 10:25 CST on Decemlier 1, 1978 at the altitude of 81 km. By 
the end of observation at lo:25 CST this layer, maintaining a constant 
scattering cross section, had descended to 75 km, leaving behind a non- 
descending layer peaking at 84 km. Note also that a weal er layer at 63 km 
maintained its altitude throughout the six-hour period. 

Individual scattering layers last from ten minutes to six hours or more 
(typically four hours) with some short-time variability, as can be seen by 
comparing the SO** and 90“* averages in Figure 5.5 or Figure 5.6. A 15-dB 
change in scattered power within five minutes is observed in Figure 5.6 around 
8:30 CST in the 80 km altitude region, a good illustration of the growth of 



Figure 5.5 Minute-by-minute scattered power profiles between 800 and 1000 
CST on April 11. 1978. 
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a whitecap. The lifetime of a whitecap is of the orJer of tens of minutes, 
suggesting that energy from internal gravity waves plays a role in generating 
the turbulence. Substantial hour-to-hour changes in the median power are 
also evident. 

Although enhancements in the scattering cross section generally occur 
in well-defined layers, there are times where these layers are absent and 
continuously enhanced scattering occurs through most of the mesosphere. One 
example of this is shown in Figure 5.8 between 10:15 and 12:15 CST on April 
13, 1978. It is likely that this broad distribution of scattered power is 
due to continuously breaking gravity waves generating turbulence over an 
extended region of the mesosphere. 

The present pulse width of the Urbana radar does not allow good 
determination of the vertical extent of most stable scattering layers, but 
the thickness is often less than the vertical resolution of approximately 
3 km. 


Estimates of the correlation time of the scattered signal have been 
found by integrating the autocorrelation function. It appears that the 
correlation time varies with height as seen in Figures 5.9 and 5.10. A 
maximum in the correlation time of about 1 sec is seen in these figures 
between 60 and 70 km. The shape of the correlation time profile (Figure 
5.9) compares well with the median half-correlation-time profile found by 
Countryman and Bowhill [1979] in the mesospheric echoes from Jicamarca, Peru. 
A maximum correlation time of approximately one second is seen in both data 
sets. The difference at the highest altitudes is due to contamination by 
the equatorial electrojet in the Jicamarca data, causing zero correlation 
time. The Urbana profile, on the other hand, levels off around 0.5 second. 

It may be significant, however, that the Jicamarca radar uses a frequency of 
50 MHz compared to the 41 MHz of the Urbana radar. 

The decrease in correlation time below 65 km, seen in both data sets, is 
probaoly due to strong photochemical relaxation of the ionization, which 
limits the lifetime of the irregularities. Below 63 km, this relaxation time 
becomes less than one second and thus affects the observed correlation time 
[Hill and Bouhill^ 1976]. Above the peak in the correlation time at about 
68 km, the decrease in correlation time with increasing altitude may be 
partly due to divergence of the radar beam. A larger scattering volume has 
less coherence in the line-of-sight velocity across the beam, thus broadening 
the spectrum and giving a shorter correlation time. It also follows that 
any broadening of the antenna pattern will shorten the correlation time at 
each altitude unless some limit is reached. 

An experiment 'as made to shed further light upon the correlation time 
variation. Three sections of the existing antenna were fed slightly out of 
phase, thus defocusing the antenna pattern. The curve marked "Wide Beam 
Antenna" in Figure 5.9 shows the height profile of the resulting correlation 
time; it has decreased markedly and has become constant with altitude. This 
suggests that the antenna beam has been broadened enough at all altitudes to 
give a scattering volume larger than the outer scale of the turbulence, thus 
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Figure 5.9 Profiles of the correlation time of the 
scattered power from the norn»al Urbana 
antenna and the modified wide beam antenna. 
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including the whole velocity spectrum. 

1116 correlation time does not seem to depend on whether a layer of 
increased scattered power is present (Figure 5.10). Particularly around 
09:30 CST is the correlation time decreasing continuously through the two 
peaks at 72 km and 81 km. lliis contrasts with the observations of Rastogi 
and Bowhill [1975], Harper and Woodman [1977] and Countryman and Bowhill [1979] 
at Jicaraarca who showed that longer correlation times are associated with 
higher scattered signal power at 70 km, where as at 77.5 km the association 
is reversed. 

Median scattered power profiles obtained between 07:00 and 17:00 CST 
are shown in Figure 5.11 for the spring (April-May 1978) and the fall (August- 
December 1978) seasons amounting to 97 and 58 hours, respectively. The peaks 
in the profiles are close to 76 km in both seasons. Most data used in this 
study were obtained between the altitude of 58.5 km and 87 km; however, some 
data from May 1978 show that very little coherent scatter is present above 
95 km. 

Figure 5.12 shows the median and 9U‘’> relative power profiles for all 
1978 data. Hie OO'c profile peaks at an altitude of 83 km compared with 76 km 
for the median. ITiis probably reflects the increased incidence of whitecaps 
at altitudes above 76 km since the 90» power level is more sensitive to these 
triur ient disturbances than is the median. It is also interesting to note 
that the discrete scattering layers noted above are not sufficiently consistent 
in altitude from day to day to appear on the cumulative profiles of Figure 5.12. 

5,7.2 Vclooities. Unlike long-period atmospheric waves, the air motions 
caused by short-period gravity waves have a substantial vertical component. 

I'hus the line-of-sight velocity observed by the Urbana radar is dominated by 
oscillations with periods of the order of ten minutes. However, on the 
scale of hours, it is known that the horizontal motions are two orders of 
magnitude larger than the vertical motions, so one may reasonably assume that 
the hourly average air motion observed is primarily due to the line-of-sight 
component of the horizontal velocity appearing in the slightly off- vertical 
antenna beam. We are therefore able to study both short-period gravity 
waves in the form of vertical velocity oscillations and long-period gravity 
waves. 

A complete study of waves in the upper atmosphere requires knowledge of 
both the vertical and horizontal wavelength as well as the period. The 
Urbana radar configuration does not allow observations of the horizontal 
wavelength so our conclusions are somewhat limited. 

A few typical examples of short-period gravity waves arc presented in 
this section. Figure 5.13 contains velocity time series for the period 8:00- 
10:00 CST on April 11, 1978. The corresponding power profiles have been shown 
in Figure 5.5. Note that the phase of the waves present throughout the period 
is a step function a*- a function of altitude, where the changes in phase 
correspond to minimums in the power profile in Figure 5.5. This is the result 
of oversampling of the signals scattered form a layered profile of turbulence. 
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Apart from this effect, it is clear that the velocity oscillations show up- 
ward phase propagation throughout all altitude regions between 8:00 CST and 
8:18 CST. This corresponds to downward propagating energy, llie vertical 
wavelength is of the order of several tens of kilometers. At 8:18 CST the 
waves in the altitude region between 78 and 84 km is replaced by nonperiodic 
motions lasting until 8:43 CST. Comparing with the power profiles in Figure 

this corresponds to the time of the large peak in the scattered power, 
interpreted as the breaking of the short-period gravity waves at a pre-existing 
scattering layer creating additional turbulence. Ibis phenomenon is quite 
frequent in the data. At lower altitudes during this time, there is a change 
from upward phase propagation to downward phase propagation involving a 
transient change of period around the 72 km level. After 08:43 CST there 
are waves with downward propagating phase throughout the whole altitude range 
from 61.5 to 84 km, representing upward propagating energy. The vertical 
wavelength is of the order of several tons of kilometers. 

Another common phenomenon is seen around 09:30 CST where the wave field 
above and below 7S km is completely unrelated. Ibis could possibly be due to 
ducting of gravity waves between a layer at 75 km and one at a higher 
altitude. 

Figuio 5,14 shows the velocities from April 13, 1978 between 1017 and 
1217 CST, corresponding to the scattered power shown in Figure 5.8, the 
scatter being fairly evenly distributed with altitude. The frequency spectrum 
of the velocity oscillations appears to vary considerably during the 2-hour 
period, perhaps due to multiple packets of gravity waves propagating through 
the region of observation. On several occasions, phase reversals with 
altitude can be seen, the best example being between 79.5 and 91 km at 
11:10 CST. This probably represents a standing-wave pattern caused by 
gravity-wave reflection at some level above 84 km. However, the short 
duration (approximately 15 minutes) of the standing-wave pattern indicates 
a short-lived reflecting layer. 

Hourly averages of the horizontal wind for the same time period as 
Figure 5.5 are shown in Figure 5.15, together with the standard deviation 
of the 1 ine-of-sight velocity and its half-correlation time. Comparing the 
horizontal velocity profiles with the corresponding power profiles in Figures 
5.5 and 5,6, we find that the peaks in the power profile occur at altitudes 
of local maxima in the S-L horizontal wind profile, or to zero vertical 
gradient in the S-li horizontal wind. 

It is, however, premature to draw conclusions about the relationship 
between wind shear and scattered power on the basis of the present data, for 
two reasons. First, nothing is known from the present data about the N-Ii 
component of the wind. Second, there may be a tendency for the wind profiles 
of Figure 5.15 to be dominated by a relatively small number of discrete 
scattering layers having some systematic relationship to the actual wind-shear 
profile, resulting in a distorted apparent wind profile. This point clearly 
needs further investigation. 

The spring and fall averages of the standard deviation of the vertical 
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Figure 5.15 Calculated hourly means of the horizontal 
velocity, the standard deviation and the 
half correlation time of the line-of-sight 
velocity profiles between 800 and 1000 CST 
on April 11, 1978. The two curves in the 
horizontal velocity profile represents the 
standard deviation profiles. 
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velocity are shown in Figure 5.16. Note the enhanced wave amplitude below 
63 km for the fall data arr* below 67.5 km for the spring data. Figure 5.17 
from Bauman and Thomae [1976], shows a local maximum in the Bnmt-Vaisala 
period around 65 km. This barrier in the Brunt-Vaisala frequency acts as a 
reflector of short-period gravity waves creating a standing-wave pattern 
below 65 km. Tlie increase in amplitude seen in Figure 5.16 aroi»nd 60 km may 
correspond to the top of this standing-wave pattern, the Brunt-Vaisala barrier 
being locat^'d at a higher altitude in the spring than in the fall. It is 
emphasized l\owever, that there are cases (such as that in Figure 5.15] where 
there is no increase below 60 Km. 

Above 60 km the wave amplitude in Figure 5.10 is essentially constant 
with altitude both for the spring and the fall data. This suggests a 
continuous absorption or reflection of wave energy with altitude, since other- 
wise the conservation of energy with altitude implies a wave amplitude that 
increases as the -0.5 power of the pressure. 

Since the Brunt-Vaisala barrier tends to suppress shorter-period waves, 
it is intuitively expected that wave spectrum above 65 km should contain 
more long-period waves than below 65 km. From Figure 5.18, it is seen that 
the correlation time increases above 63 km in the spring and above 60 km in 
the fall probably corresponding to the lower edge of the Brunt-Vaisala barrier. 
Below 60 km the correlation time indicates a wave period of about 4.5 to 5 
minutes, close to t‘'e Brunt-Vaisala period. Between 66 and 80 km the average 
wave period is constant at around 8 minutes in the spring and about 12 minutes 
in the fall. The reason for this difference in wave periods is unclear, 
but could be due to variations in the strength of the Brunt-Vaisala barrier. 
Above about 80 km the correlation time starts falling off, indicating a 
cascading of energy from long-period gravity waves in to short-period gravity 
waves and turbulence above 80 km. This interpretation is consistent with the 
data of figure 5.13 showing breaking waves or whitecaps above 78 km. At 
87 km both the spring and the fall correlation times again approach a value 
corresponding to the local Brunt-Vaisala period of about 5 minutes. 

5.7.3 Disauesion and conclusion. Scattered power in the mesosphere 
tends to occur in persistent layers with a separation of a few km, as seen 
in Figure 5.5. It is tempting to assume that this separation is related to 
the vertical wavelength of tides or long-period gravity wav.-js. However, the 
nearly constant altitude of most of these scattering layers then requires a 
standing-wave pattern, caused by wave reflection at an altitude above 90 km. 
Reflection of tidal components has been reported by several groups [Fallous 
ct at., 1975; Tsuda et at., 1979]. Vertical wavelengths as deduced from the 
power profiles range from 12 to 36 km; or 6 or 18 km if there is only one 
layer per wavelength. 

Scattering profiles such as Figure 5.7 can be interpreted as a combination 
of a standing and a descending tidal component resulting from partial reflec- 
tion of a tidal wave. 

Since turbulence is related to velocity gradients we expect tides with 
short vertical wavelengths to be the dominant cause, since long vertical wave- 
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Figure 5.18 Seasonal correlation time profiles for the 
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lengths give very small gradients even if the amplitude is large. Both 
b'elloua et at. [1975] and Tauda et al. [1979] have shown that observed 
diurnal tide has much shorter vertical wavelength than the observed semi- 
diurnal tide. Arsuming that the descending tide in Figure 5.7 is diurnal, 
we have a vertical wavelength of 24 km which corresponds well with the 
expected vertical wavelength of the 5]! diurnal mode at the altitude of 80 km. 

The standing-wave pattern in Figure 5.6 is most likely due to a diumal 
or semidiurnal tide with vertical wavelength of 18 km. This compares well 
with the south-east component of the horizontal wind in Figure 5.15 which 
shows a dominant vertical wavelength of 18 km. 

Although stable layers usually dominate the scattered-power profile, 
there are still substantial variations in amplitude and scattering height 
(Figure 5.5). Since the short-term variations in scattered power have 
periods (about 10 min) comparable to those of the velocity fluctuations, 
we believe tha^ these short-term varia lions are caused by gravity wave 
modifications ci the wind profile, thus varying the amplitude and altitude 
of maximum shear, and hence the turbulence. However, we have been unable 
as yet to establish a detailed relationship between wave velocity and 
scattered power. 

Tlie amounf of scattered power depends also on the background electron- 
density gradient, which is cdvecteJ by the turbulent velocity field, subject 
to the effects • f relaxation {Hill and Bowhillt 1975]. Tlie mutual interaction 
between these effects is quite complex, and the electron-density gradient 
IS quite variable from day to day. However, the averaged electron-density- 
gradient profile (Figure 5.19, a component from Meahtly et al. [1972]) 
increases with altitude except for local maxima at the bottom ledge of the 
D and E layers respectively. In addition, the photochemical relaxation 
rate is higher at lower altitudes, thus favoring gr^'ater scattering from 
higher altitudes. 

Below 70 km the decrease in turbulent energy with increasing altitude, 
at the wave number of the Urbana radar, in linear [Raatogi and Bowhillt 
1976b]. However, the energy-spectrum of the turbulence falls off exponentially 
as a function of altitude above approximately 70 km. Furthermore, above 80 
km the (^-folding distance of the exponential decay decreases, indicating an 
even greater decrease of turbulence with altitude. Thus, we believe the 
decrease in observed scattered power with altitude to be due to the sharp 
decrease in turbulent energy overriding the increase in electron-density 
gradient above about 70 km. 

We have seen that vertical wavelength of short-period gravity waves is 
large, 30 km or more, for the gravivy waves with periods around 10 min. and 
that the averaged vertical wave-velocity is more or less constant above 65 
km. This agrees well with computation by Bowman and Thomas [1976] who 
found that the vertical gravity wave motion reaches a maximum around 60 km 
whereas tlie horizontal wave motion has maximum at 85 km. Thus, we conclude 
that the observations of the gravity wave profile confirms the existence 
of the Bru/it-Vaisala barrier between 60 and 65 km. 



Figure 5.19 Profiles of mean electron density and electron density gradient 
in the mesosphere. 
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In sununary, the data suggests that the VliF scattered power profile from 
the mesosphere over Urbana depends on the existence of short-vortical-wave- 
length tidal components to generate turbulence. In addition, short-period 
altitude and aji.plitude variations are believed to be due to gravity wave 
interacting with the tidal components. These short-period gravity-wave 
profiles show a variation that is consistent with theoretical predictions of 
Bauman and Thomas [1976]. 
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6. LASER RADAR 

A laser-radar syst m has been designed and constructed to measure the 
altitude distribution of atmospheric sodium in the 80- to 100-km altitude 
region. Ibc system consists of a high-power pulsed dye laser, large aperture 
receiving telescope, and photon- counting and signal-processing equipment. 
Coordination of the project has been provided by Prof. C. F. Sechrist, Jr. 

Also involved are Prof. C. S. Gardner who has been concerned with the receiving 
system development, and Prof. H. Merkelo of the Quantum Electronics Research 
Laboratory who has offered assistance in the construction and development of 
the dye laser. This research is supported in part by the National Science 
Foundation under grant ATM 79-20726, and in part by the National Aeronautics 
and Space Administration under grant NGR 14-005-181 

6.1 Tuo- Dimensional Signal Pvaceeuing 

Theoretically, the laser pulse width, receiver range gate width, data 
integration time and sampling interval place the fundamental limits on the 
spatial and temporal resolution of lidar profiles. However, in practice 
the resolution of sodium lidar data is usually limited by shot noise. Spatial 
filtering has been shown to be an effective technique for enhancing the 
spatial resolution of sodium lidar profiles by reducing shot noise in the 
photoc.ount data \Rox:>lett et al.^ 1978; lioiilett and Gardner, 1979]. During 
this past year J. 1). Shelton and C. S. Gardner have developed a two- 
dimensional signal processing technique which utilizes spatial and temporal 
filtering to further reduce shot noise and increase profile resolution 
[Gardner and Shelton^ 1981; MihavioBt 1979]. The 2-D technique has been 
applied to our sodium lidar data since the fall of 1979 and has revealed 
some interesting small-scale short-term variations in the layer. 

The most critical aspect of the 2-D filtering process is the determin- 
ation of the appropriate spatial and temporal cutoff frequencies. This is 
accomplished by examining the spatial and tempore 1 periodograms of the photo- 
count data. First, the average background noise level is subtracted from the 
raw photocounts comprising each sodium profile. The spatial periodogram of 
each profile is computed using a fast Fourier transform (FFT) algorithm. Then 
the periodograms of all the profiles collected during the night are averaged. 
The average spatial periodogram for the 71 profiles collected on the night 
of October 28-29, 1979 is plotted in Figure 6.1a. The dashed line is the 
average shot noise level which is approximately 3t dB below the spectral peak 
at wave number zero. The dashed curves are one standard deviation above and 
below the average periodogram. Because of the large number of profiles used 
to compute the average, the error is quite small. Notice that the average 
periodogram falls into the shot noise at a wave number near 0.3 km"i. 

The temporal periodogram is more difficult to comyjute and analyze for 
two reasons. First, temooral variations of the si^dium layer at a fixed 
altitude are usually a very small percentage of the average sodium background. 
The second difficulty arises from the shaip discontinuities in the temporal 
data at the beginning and end of the observation period. The sharp 
discontinuities along with the high average level dominate the temporal 
periodogram and overshadow spectral componenets due to actual temporal varia- 
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Iigurc 6.1 Spatial and temporal poriodograms of the sodium lidar 
photocount data collected during the night of October 
28-29, 1980. 

(a) The cutoff frequency for the temporally filtered 
data is 0.05 min' * . 

(b) The cutoff wave number for the spatially filte.cd 
data is 0.44 km' ^ . 
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tions in the sodiiiin layer. This problem can be ameliorated by subtracting 
the average sodi'im background level from each range bin and then multiplying 
the resulting tein;>oral data by a window function to eliminate the discon- 
tinuities at the beginning and end of the observation period. The temporal 
periodogram for each range bln can then be computed by taking the FFT of the 
windowed data. 

Unfortunately, the periodogram provides a very noisy estimate of the 
signal spectrum. The noise can be reduced and a more reliable estimate of 
the spectral features present in the temporal spectium can be obtained by 
averaging the periodograms for all range bins in the sodium layer. The 
dominant spectral <^catures present throughout the layer during the observation 
period will be retained in the average temporal periodogram. Because we 
are computing the temporal periodogram for the purpose of choosing an 
appropriate filter cutoff frequency and arc not interested in the detailed 
spectral structure, the slight distortion introduced by windowing the data 
and averaging the periodograms is nut significant. 

Tlie average temporal periodogram for the 71 profiles collected during 
the night of October 28-29, 1979 is plotted in Figure 6.1b. The peak 
signal level near temporal frequency 0.003 min"^ is approximately 10 dB 
above the shot noise level indicated by the dashed line. Since the signal 
falls into the shot noise near the frequency 0.03 min*^, a lowpass filter 
with a cutoff of 0.05 min‘^ will be used to process the data. Because the 
temporal filter is applied directly to the raw photocount t’ota, a finite 
impulse response filter is required to prevent distortion which could be 
introduced by the short discontinuities in the data at the beginning and 
end of the observation period. We have chosen to use a lowpass filter with 
an impulse response which is identical to the Hamming window function. 

Temporal filtering can result in a significant reduction in shot noise. 

To illustrate, the October 28-29 data were lowpass filtered with a cutoff 
frequency of 0.03 min"^. Here we define the cutoff frequency as the first 
cero in the OFT of the Hamming window function. Tlie average spatial period- 
ogram of the temporally filtered data is plotted in Figure 6.1a. By comparing 
the spatial periodogram of the temporally filtered data with the periodogram 
of the raw photocount data also plotted in Figure 6.1a, it is clearly seen 
that temporal filtering has reduced the shot noise level by more thtui 10 dB. 
Spectral features which were completely hidden by the shot noise in the raw 
data are readily observed in temporaily filtered data. In a similar manner, 
spatial filtering can reduce the shot noise level in the temporal periodogram. 
This is illustrated in Figure 6.1b where the average temporal periodograms are 
plotted for the raw data an(' spatially filtered data collected on October 28- 
29. For this particular data set and choices for the spatial and temporal 
cutoff frequencies, the approximately 3 dB reduction in the shot noise level 
due to spatial filtering is not as dramatic as the 10 dB reduction due to 
temporal filtering. 

Afver the appropriate cutoff frequencies are chosen by examining the 
spatial and Lcmporal periodograms, the raw sodium profiles are two-dimensional ly 
filtered. The filtered data are then corrected for the inverse square altitude 
dependence of the lidar signal and the results plotted. In Figure 6.2 are 
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Figure 6.2 Sodium density profiles obtained by 2-0 filtering the 
photocount datd collected on October 28-29, 1980. The 
temporal cutoff frequency is 0.03 min"^ and the spatial 
cutoff wave number is 0.44 km" ^ ; 2-1) plot with the 
horizontal axis represent, ng relative sodium concentration 
and time. 
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plotted the 2-D processed sodium profiles for the data collected on October 
28-2P, The time inf'rer.cnt between profiles is 15 minutes. Figure 6.3 is a 
3-D plot of this same data. These plots clearly show the propagation of 
Wi. e-like features through he layer. In the 3-D plot the phase progression 
O’’ vave is quit<5 apparent on the bottom side of the layer. Notice also 

the 3-D plot that the sodium concentration near the layer peak at 90 k.n 
itude is undulating. The leabon for this undulation is evident in 
igurc 6.2. The layer peak mrvos fix)m an altitude of 91 km at 2233 CST to 
slightly below 90 Vm a^ 0048 CST aid then returns to 91 km at 0233 CST. 

Rowlett and Gardner [1979], Mihaviae [1979], and Gardner and Shelton 
[1981] provide more detailed descriptions of the spatial and 2-D filtering 
techniques . 

6.2 Computer-Model Studies of the Sodium Density Response to Gravity-Wave 

Perturbations 

Tfie density response of the sodium layer to gravity waves is being 
investigated theoretically and experimentally by J. D. Shelton. In general, 
the response of an atmospheric layer to a gravity wave may be determined by 
examin ng the equation of continuity. This differential equation indicates 
that the velocity perturbations associated with a gravity wave will induce 
density perturbations in the sodium layer. Assuming the diffusion time for 
sodium is greater than the gravity wave period being considered, the sodium 
velocity perturbations will be equal to the velocity perturbations in the 
background atmosphere caused by the passage of this wave These velocity 
perturbations may be expressed in terras of atmospheric density variations 
by utilizing the polarization relations for gravity waves. As a result, 
the solution of the continuity equation will express the density response of 
the sodium layer in terms of atmospheric density perturbations caused by 
gravity waves. A perturbation series solution of the continuity equation 
may be found. This solution consists of a sum of terms. The first term in 
this solution describes the linear layer response while the nonlinear portion 
of the layer response is contained in the higher order terms. In general, 
the layer response is dominated by the linear term. 

Several features are associated with the linear term. These include: 

(a) larger density perturbations below the layer peak than above it. 

(b) density perturbations with the same temporal frequency and hori- 
zontal wave number as the gravity wave. 

(c) a point of phase reversal in the linear response occurring slightly 
above the layer peak. 

Near the point of phase reversal, the lir. ar response is small and higher 
order terms may make significant con .ributions to the total layer response. 

The frequencies of these higher order terms are harmonics of the 
gravity wave frequency. For example, perturbations associated with the second 
ord.ar term have twice the frequency of the gravity wave. In addition, if 
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Sodium density profiles obtained by 2-D 
photocount data collected on October 28 29, . 

temporal cutoff frequency is 0.03 min and the spati 

clZff wave number is 0.44 km' ^ 3-D plot of the same 
profiles plotted in Figure 6.2. 
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the perturbations induced in atmospheric density by the gravity wave are 
small, the magnitude of these high order terms rapidly decreases with 
increasing order. Figure 6.4 contains a plot of the envelopes of simulated 
density perturbations associated with the first and second terms of the 
perturiation series solution and their sum. A Gaussian layer and 5« peak- 
to-pcak variations in atmospheric density have been assumed. It is evident 
that a double frequency component in the sodium density perturbations should 
become dominant in the 90.3 to 91.2 km altitude range. This type of behavior 
has been observed in experimental data. Figure 6.5a contains a plot of the 
envelope of density variations as well as the average sodium density profile 
observed on October 28-29, 1979 Figure 6.5 contains a corresponding plot 
of simulated data assuming the odium density profile in the absence of gravity 
waves is the same as the average density profile given in Figure 6.5a. It is 
seen that good correlation exists between the curves of maximum and minimum 
density in experimental and simulated data. 

6.3 Neu Laser for Lodiu/'i LI PAR 


A new fl ash ljunp- pumped dye laser for use as the transmitter in the Urbana 
LIDAR system has recently been obtained. The required specifications for this 
laser arc given below: 


Center frequency 
Linewidth, FWIW 
energy per pulse 


5890A 

O.OIA 


SOmJ 


Config. B 
5890A 
0.002A 
20mJ 


Configuration B is identical to Configuration A with the addition of an 
electronically tunable etalon to further reduce the linewidth. To have the 
capability to measure the linewidth of the new laser a spherical mirror 
scanning interferometer was designed and constructed by R. Nilsson. Testing 
was done by inserting a "chopper" at the output of a thermally stabilized 
MeNc alignment laser to simulate a pulsed laser, ilie purpose of the testing 
was to determine that the linewidth measurement system performed adequately 
;md to establish that interferometer instabilities were not great enough to 
interfere with the linewidth measurement. There is little broadening of 
measured linewidth due to any frequency instabilities. Tliis implies that 
if this interferometer is used to measure the linewidth of the new laser 
there will be negligible broadening due to instabilities in the interferometer. 


When attempts were made to use this interferometer to measure the line- 
width of the new laser it provided to be so difficult to align and the out- 
put fringes of such poor quality that meaningful results were not obtained 
with it. The spherical mirrors were replaced with a pair of flat 85"> mirrors. 
Iliese mirrors had a measured finesse of 8.8 and, {"or the chosen interferometer 
spacing this yielded a value of Av„,in output obtained using 
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Figurc(v-1 Predicted envelope of normalized density perturbations due 
to first three terms of perturbation series solution. 5“ 
peak-to-peak background density variations are assumed. 
Gaussian layer described by = N exp [-('’-90]' /IS] is 

assumed. 
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Figure ti.5 (a) Fxperimental ly observed envelope of density 
variations on October 28-29, 1979. Solid line is 
average density profile. (b1 Predicted density 
variations due to first- and second-order terms. 
Solid line is the assumed steady-state density 
prof i le. 
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those mirrors made possible a linewidth measurement, but with a value of 
AV|^i_ of yso Wiz, the convolution of the interferometer frequency trans- 
mission function and the laser linewidth cannot be neglected; it could have 
been neglected had the high finesse of the spherical mirrors been available. 
Knowing the value of the laser linewidth can be calculated from the 

measured linewidth. Using this knowledge a linewidth of 900 Miz was measured 
at an energy of .SO mJ/pulse for Configuration A and a linewidth of 400 Miz 
(O.OOSA) at an energy of 20 nU/pulse for Configuration B. 

From the measured values of the laser linewidth it is concluded that 
Configuration A met all of the required specifications but that Configuration 
B did not. To be sure that the measured value of the laser linewidth was 
accurate, considering the difficulties encountered with the interferometer, 
the linewidth measurement was repeated with a Spectra- Physics Model 450-03 
scanning interferometer, which is now available. The Spectra-Physics inter- 
ferometer is a completely sealed unit. Mirrors were aligned at the factory 
and a photodetector is included. With a finesse of 200 and an FSR of 2 GHz, 
a bandwidth of 10 MHz results; this is ^imall enough to neglect the convolution 
effect of the laser linewidth with the interferometer bandpass. The repeated 
linewidth measurement with the Spectra-Physics interferometer again resulted 
in a bandwidth of 400 MIz (0.005 ?l) for Configuration B at the specified 
energy of 20 mJ/pulse. 

At this time further work is being done by the laser manufacturer. 

Candela Corji., in an attem^jt to r.cke the laser conform to the specifications 
of Configuration B. 


. 4 LSI ‘21 Mi ovoootripu ter 

ITie University of Illinois sodium lidar system presently uses a 2-80 
microcomputer to handle its receiver operations and data collection. Soon, 
an LSI-11 microcomputer will replace the 2-80 for several reasons which are 
outlined below. 

A pulsed dye laser (manufactured by Candela Corporation; see Table 6.1 
for laser parameters) is tuned to the D 2 sodium resonance line at a wavelength 
of 589 nm. Refer to the block diagram of Figure 6.6. This is done by 
directing a small portion of the laser output to a sodium vapor cell which 
serves as a wavelength monitor. The sodium vapor will fluoresce when a 
properly tuned laser beam passes through the cell. Wlien the lidar system is 
properly aligned, the receiving telescope collects a fraction of the laser 
energy which has been backscattcred from the sodium layer. The telescope 
filters the collected light to minimize background noise and then focuses the 
filtered signal onto the photocathode of the photomultiplier tube (PMT) . The 
PMT is thermoelectrical ly cooled to -20“C in order to reduce dark nol.-^e. Ilie 
discriminator generates logic-level pulses whenever a photon is detected by 
the PMT. The discriminator is connected to the photon-counter processor which 
contains a high-speed counter and semiconductor memory. 

The microcomputer interfaces wit.'i the lidar receiving system directly 
through the photon-counter processor. After the laser has fired, the computer 
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SODIUM LAYER (80-90 Km ) 



figure 6.6 Block diagram of Urbana sodium lidar system. 
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instructs the processor to sample the photon> counter periodically and store 
samples in a semiconductor memory. 256 samples are obtained for each I'iser 
shot. ITie sampling interval can be varied from 1 to 256 us. The one usee 
sampling interval corresponds to a 150-m range resolution. When all 256 
samples have been collected, the data are transferred from the processor 
memory to ihe computer. The computer subtracts successive samples to obtain 
the number of photons counted within each sampling interval or range bin. 

The photon counter continually increments and is reset only when it overflows. 
The photocounts are accumulated for each range bin until enough photons have 
been .wcted to construct a sodium profile. Typically 50 to 200 laser shots 
are required to obtain a profile. The raw data profiles are processed off- 
line or a CYBER- 175 computer to obtain profiles of the sodiuic density versus 
altitude and time. 

Presently a Z-80 performs the tasks of the microcomputer described above 
in addition to more minor tasks such as: (1) setting time duration between 
laser shots (typically 2 to 3 sec.), (2) aiding laser and receiver beam 
alignment, and (3) storing raw sodium profiles on cassette tape. 

The LSl-11 microcomputer was chosen to replace the Z-80 because of the 
much faster processor ajid larger memory cap.*icity. The faster processor 
should permit collection of lidar data at the maximum repetition rate of 
the new Candela laser (which is iw Hz or once every 100 msec). The Z-80 will 
permit us to collect data at only a 2 Hz rate or slower. The faster processor 
coupled with a larger memory (20K vs. the Z-80's 8K) will also allow some 
basic on-line signal processing, thus permitting real-time observation of the 
sodium layer profile. Also, the LSI-11 will use a dual floppy disk-drive 
system to load programs and store raw profiles. The floppy disks provide a 
much faster means of transferring data to the CYBER 175 computer. Additional 
LSl-11 tasks will include monitoring the laser parameters of wavelength, line- 
width, and energy per pulse as well as controlling the laser tuning and 
firing. 

Die complete LSI-il package is centered around the OEC LSl-11/2 processor 
and housed in an MDB Systems backplane operating unit; the system also includes 
two DSD 210 floppy disk-drive units, 20K words of memory and two peripherals, 
one DEC DLVll-J, 4-port asynchronous serial line unit and one DEC DRVll parallel 
line unit. 

6.5 Reaeiving System Improvements 

A thermoelectrical ly cooled housing has been acquired for cooling the 
receiver photomultiplier tube. Fhis unit, a Products For Research model 
TE-206 TS-RF housing, consists of an air h*^at exchanged housing, associated 
power supply, and a double pane pyrex window. An integral part of this window 
is a 1.5-inch focal length lens which focuses collimated light from the 
receiver optics onto the PMT photocathode. This unit replaces a unit which 
required liquid nitrogen vapors for cooling, thus alleviating the need for 
an reservior of liquid nitrogen. Both housings accommodate the RCA C31034A 
photomultiplier tube. 
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6.(1 l)lavelK.riijth Meaeurcfncnt 

The feasibility of utilizing a Fizoau waveraeter to accurately monitor 
the wavelength of the Urbana LIDAR transmitter output beam was investigated 
by 0. Burlet. Tliis instrument is based on a Fizeau or "optical wedge" 
interferometer, the fringe pattern of which is digitized and stored in a 
small computer. TJiis information is then processed to detemine the wave- 
length of the interfering light. 

The lizeau wavemeter divides the collimated input beam into two overlapping 
beams of equal intensity, propagating not quite parallel to each other. TTie 
i'*sulting uniformly space Fizeau fringe pattern is then sjunpled periodically 
by an array of photodiodes which digitizes the pattern and sends it into the 
computer memory for processing. 

i'rocessing consists of the following; the computer determines the spatial 
period of the patteiTi, and the location of a fringe intensity niinimiaii, Tlie 
interferometer wedge angle whidi is a known constant then allows the spatial 
patteni to be converted into a first approximation of the wavelength. 'Hie 
uncertainty here must be less than the free spectral range of the interferometer. 

The computer then determines the interferometer path difference at the 
previously determined fringe minimum, utilizing the wedge angle and wedge 
spacing at a known reference position. The path difference at a fringe 
minimum contains an integer number of wavelengths. The computer, therefore, 
calculates the ratio of the path difference to the approximate wavelengt’.i 
determined abrve, which gives the interferometer order number plus a 
remainder. Hie remainder is then used to calculrce the ratio of the remainder 
to the order niimber, this then being applied to the approximate wavelength, 
yielding a new estimate of the wavelength. This procouure is repeated until 
the solution converges to a "best value" of the wavelength limited only by 
the resolution of the measurement of the path difference at fringe minimum. 

The study concluded the wavemeter is a feasible method of monitoring 
the output provided an enclosed and evacuated chamber for the interferometer 
can be easily constructed to eliminate dispersion and scattering of the 
interfering becuns (see Figure 6.7J. 
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[hiring the regular Fall and Spring aaMatera, a aeainar series in 
aeronomy is offered as graduate course EE 490. The fallowing seninars were 
conducted during the period covered by this report under the supervision of 
Prof. C. F. Sechrist, Jr. and Prof. S. K. Avery. 


October 7 

Mr. Carl Stehle 
Graduate Assistant 
Department of Computer Science 
University of Illinois 

A Neutral Model of the 
Thermosphere in Magnetic 
Coordinates 

October 14 

Dr. Jane Fox 

Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

A Stratospheric Chemical 
Instability 

October 17* 

Dr. J. Bufton & Mr. J. Abshire 
NASA'Goddard Space Flight Center 
Greenbelt, F® 

Laser Radar and Ranging 
Projects at NASA- Goddard 
Space Flight Center 

October 21 

Dr. Wayne M. Wendland 
Illinois State Water Survey 
University of Illinois 

Advection Patterns In- 
ferred from Mt. St. Helens 
Volcanic Plume 

Noven4)er 4 

Mr. Gary A. Jones 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Planetary Waves and Nitric 
Oxide Transport 

November 11 

Dr. P. M. Dolas 

Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Gravity Waves in the 
Equatorial Mesosphere 

November 25 

Prof. S. K. Avery 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Correlative Measurements 
Using *he Urbana lidar 
and Meteor Radar 

December 2 

Prot. S. A. Bowhill 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Radar Signal Processing by 
Microcomputer: I. Tlie 

FORTH Language 
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Prof. S. A. Bowhill 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Radar Signal Processing by 
Microcomputer: II. Imple- 

mentation of a Coherent 
Scatter System 
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Laser Applications: The Early 
Days of Uranium Isotope 
Separation 


February 3** Albert G. Brejcha The Voyager Program 

Jet Propulsion Laboratory 
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Prof. S. A. Bowhi'l 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 

Prof. Gemot Metze 
Dept, of Electrical Engineering 
Aeronomy Laboratory 
University of Illinois 
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